Effective crystal field and Fermi surface topology: a comparison of d- and dp-orbital 
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The effect of electronic correlations to enhance or reduce the effective crystal field in multi-orbital 
correlated materials can be crucial in determining the topology of the Fermi surface and, hence, the 
physical properties of these systems. In this respect, recent local density approximation (LDA) plus 
dynamical mean-field theory (DMFT) studies of Ni-based heterostructure have shown contradicting 
results, depending on whether the less correlated p orbitals are included or not. We clarify the origin 
of this problem and identify the key parameters controlling the Fermi surface properties of these 
systems. A particularly important one is the filling of the identified d-orbitals: in the dp-calculation 
this is typically larger so that Hund's exchange leads to a different local magnetic moment for the 
dp-model. 

PACS numbers: 71.10.-w, 71.27.+a, 73.40.-c 
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I. INTRODUCTION 

Correlated electronic systems display some of the most 
fascinating phenomena in the solid state physics. One of 
their typical characteristic is a strong sensitivity to small 
changes of external control parameters. Hence, a precise 
understanding of the underlying physics and of the piv- 
otal parameters controlling the observed phenomenology 
represents a crucial goal in contemporary condensed mat- 
ter research, also in light of possible applications beyond 
the purely scientific context. 

The intrinsic complexity of many-body physics pre- 
vents an exact ab-initio theoretical description of corre- 
lated materials. In fact, even one of the most basic mod- 
els for electronic correlations, i.e. the Hubbard modelfll| 
where only the local part of the Coulomb interaction is 
retained, cannot be exactly solved in the relevant cases of 
two or three dimensions. However, a great step forward 
in the theoretical analysis of electronic correlation in con- 
densed matter was achieved in the last decades by means 
of dynamical mean- field theory (DMFT)[i, [|, and, for 
realistic material calculations, by its merger with ab- 
initio density functional approaches (LDA+DMFT) [i] . 

From the theoretical point of view, DMFT-based 
methods can be viewed as a quantum extension of the 
classical mean-field approaches. Hence, the application 
of DMFT implies neglecting all non-local spatial corre- 
lations beyond the standard (static) mean-field theory 
level. At the same time, DMFT as well as LDA-l-DMFT, 
allows for a very accurate (and non-perturbative) treat- 
ment of the most relevant part of the electronic correla- 
tions stemming from the (single or multi-orbital) Hub- 
bard interaction, i.e. their purely local part. The ability 
of DMFT to capture these local quantum fluctuations is 
one of the keys behind its success in addressing many 
open questions in the physics of strongly correlated ma- 
terials. Among those jWe recall the pioneering DMFT de- 
scription of the MottfSI] metal- insulator transition (MIT) 



in ¥203(1, @], of the S phase of PuJ,], of the correlation 
effectsQ in Fe and Ni, of the volume collapse in Ce[l], of 
the unconventional pairing mechanism of superconduc- 
tivity in fullerenes[l(j|: most recently DMFT has been 
also successfully applied to the analysis of the occurrence 
of kinks in the self-energy jllj and in the specific heat[T2j 
of particular vanadates, such as SrVOa, and LiV204, as 
well as of the spectral and magnetic properties'^ of Fe- 
based superconductors. 

Furthermore, it should also be recalled here, that 
DMFT is a very flexible scheme, whose application is 
possible also beyond the standard case of bulk corre- 
lated systems. In fact, DMFT-based methods have been 
recently used to study correlated nanoscopic fl^ and 
hetero-structures[15- 17]. For the latter case, we want to 
focus here, in particular, on the theoretical predictions 
for the Fermi-surface properties of layered Ni-based het- 
erostructures. The application of LDA-l-DMFT to this 
problem, and more specifically, to the case of a 1 : 1 lay- 
ered LaNiOa/LaAlOa heterostructure has raised a con- 
siderable interest, as the DMFT results of Refs. [H El] 
clearly prospect the possibility to drive the Fermi sur- 
face "topology" of these materials very close to the one 
of the high-temperature superconducting cuprates. In 
fact, the electronic structure in the bulk Nickelates, such 
as Ri_j;Sra;Ni04, is typically characterized by two bands 
crossing the Fermi level [l^, which arise from the two eg 



orbitals of Ni (the "planar" x 



and "axial" 82 



orbital). However, by growing heterostructures with 
planes of LaNiOa intercalated with insulating planes of 
LaAlOs and on substrates providing an epitaxial strain, 
such as SrTiOs or PrScOa, the energetic configuration 



of the 



will be correspondingly disfavored 20] , 



which corresponds to a (positive) crystal field splitting 
^CF — Esz^-r^ — £x^~y^ > among the two Cg-Ni or- 
bitals. In this situation, LDA -f DMFT calculations have 
shown that the inclusion of the correlation effects will 
always increase the original (LDA) crystal field splitting 
AcF among the two eg orbitals, leading, eventually, to a 
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significant change of tlie "topology" of the Fermi surface, 
i.e. to a situation in which only one band, with predomi- 
nant — y^ character crosses the Fermi level. Hence, ac- 
cording to these LDA-t-DMFT calculations, in Ni-based 
heterostructurcs it would be possible to "artificially" re- 
alize the electronic configuration of the high-temperature 
superconducting cuprates (i.e., single, almost half-filled 
orbital with — symmetry close the Fermi level), e.g. 
by modulating the strain through changing the substrate 
or the insulating layer. As the control of the low-Fermi 
surface properties represents an essential ingredient for 
novel, alternative, realizations of high-temperature su- 
perconductivity, the importance of having highly accu- 
rate LDA-I-DMFT predictions becomes a crucial factor 
for engineering new materials. 

In contrast to these impressive applications of 
DMFT-based methods, one important aspect should be 
stressed here: most of the above-mentioned DMFT or 
LDA-I-DMFT calculations have been performed by ap- 
plying DMFT to a "reduced" basis-set including only the 
most important correlated c?-orbitals close to the Fermi 
level. 

In the last years, several LDA-I-DMFT calculations 
included also the less correlated, typically ligand p, or- 
bitals. In most of the cases, these additional [p) orbitals 
were way more extended than the correlated ones (e.g. 
d), and, hence, the local Coulomb interaction between 
electrons occupying these additional orbitals (C/pp) and 
between electrons on different orbital manifolds {Upd) 
was either completely neglected, or, at most, treated 
at the Hartree level [2 1|. Despite this approximation, 
it is quite reasonable to expect that the accuracy of 
LDA-I-DMFT calculations performed in an enlarged (say: 
dp) basis-set to be higher than the corresponding ones in 
the restricted d manifold. In fact, quite generally: (i) 
performing a local (Wannier[22|, NMT0[2ii], etc.) pro- 
jection on an enlarged basis set allows for a better local- 
ization of the orbitals of the correlated manifold (as the 
dp-hopping processes are now explicitly included in the 
model); and (ii) the possibility of describing explicitly 
charge-transfer processes between the d and p orbitals 
makes evidently closer the theoretical modeling to the 
actual material physics. 

Notwithstanding these quite general arguments, when 
actually performing LDA+DMFT calculations on en- 
larged dp basis-sets the improvement w.r.t. older 
LDA-I-DMFT calculations restricted to purely d-basis is 
not always evident. Comparing with experiments, there 
are cases in which the use of the larger dp basis-set makes 
the accuracy of the theoretical description surprisingly 
worse. Without attempting to give a complete review 
here, we recall that LDA-I-DMFT calculations including 
p orbitals have improved the descriptions of the insu- 
lating behavior of NiOjH] and of the MIT in NiSsfll 
w.r.t. d only calculations f26j. Also quite accurate re- 
sults have been obtained for one- and two-particle prop- 
erties of cobaltates (such as SrCoOaf??! and LaCoOs 
[28j ) and for the well-known class of iron-pnictides and 



calchogenides [1^ . 

In contrast to the aforementioned successful applica- 
tions of dp calculations, in other, equally important, cases 
the dp LDA-I-DMFT results are in partial or total con- 
tradiction with the d-only calculations, and/or with the 
experimental findings: No MIT in V2O3 was found up to 
very large values of the Coulomb interaction if the oxygen 
p orbitals are included in LDA-I-DMFT calculations 30] . 
More recently also the Mott-Hubbard insulating phase 
of La2Cu04 and LaNiOa was reported to be missing in 
the dp framework [sif . while these materials are found 
to be insulating in d-only calculations for plausible val- 
ues of the dd interaction. These discrepancies between 
d-only and dp-calculations regarding the Mott-Hubbard 
MIT have already raised a discussion in the recent lit- 
erature. In Ref. [31[, non-local correlations neglected by 
DMFT have been considered as a cause of the discrep- 
ancy. In fact, a major role in the determining the onset 
of insulating cases is likely played also by non-local cor- 
relations, especially in the two-dimensional cuprates. 

However, there are also other discrepancies, whose 
discussion will be at the center of this work, which 
can be hardly attributed to the effects of non local- 
correlations. These discrepancies are observed for sys- 
tems of more than one correlated d orbital and in broad 
parameter regimes (including high-T), where effects be- 
yond DMFT should not play any crucial role. In par- 
ticular, a striking disagreement between d-only and dp- 
calculations was reported for the above-mentioned case 
of Ni-based heterostructures. In fact, LDA-I-DMFT cal- 
culations performed including also the p-Oxygen orbitals 
have shownjs^] exactly the opposite trend w.r.t. the pre- 
vious ones: Even in presence of a favorable crystal field 
splitting A(7i? — 63^2 — £a;2_j,2 > at the starting 
(LDA) level, the net effect of the Hubbard interaction 
was always to reduce the orbital polarization by filling 
back the — orbital, which would prohibit "de facto" 
any possibility of realizing the cuprate conditions for the 
onset of an unconventional superconductivity. 

This second kind of discrepancies between 
LDA+DMFT performed with different (d-only vs. 
dp) basis-sets well illustrated by the contradicting 
results for the Ni-based heterostructures, raises a quite 
general question about the proper use and interpretation 
of the growing number of LDA-I-DMFT calculations on 
extended basis sets. The correct determination of the 
orbital polarization and Fermi surface properties is of 
great importance for future calculations of increasingly 
complex materials. 

In this paper we aim at understanding the relations 
between the results of LDA-f DMFT calculations on dif- 
ferent basis-sets for the Fermi surface properties of multi- 
orbital systems, and, ultimately, the origin of the qualita- 
tive discrepancies, observed by following the standard im- 
plementation of the algorithms in different basis-sets. For 
this purpose, it is of primary importance to disentangle 
the main, qualitative, trends from the specific features of 
a selected case. Hence, we will select two (corresponding) 
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minimal models for d-only and dp calculations. They are 
intended to capture the above-mentioned discrepancy for 
the prediction of the Fermi surfaces of correlated multi- 
orbital systems, and allow for a systematic study of this 
problem. 

The scheme of the paper is the following: In Sec. 
2 we introduce the simplified models for capturing the 
discrepancy between the theoretical predictions for the 
Fermi surface properties of multi-orbital d and dp sys- 
tems, such as those reported for realistic calculations of 
Ni-based materials. In Sec. 3, we analyze the effective 
crystal field splitting and discrepancies between d- and 
dp-model at for a fixed set of interaction parameters. At 
the end of Sec. 3, we also provide an analysis of the role of 
Hund's exchange J. In Sec. 4 we study the origin of the 
observed discrepancies by analyzing the dependence on 
the d-electron-density, and we discuss its possible relation 
with the crossoverjSSj from high-spin (Hund's regime) to 
low-spin (CF regime). Finally, in Sec. 5 we summarize 
our results. 



II. MODELS AND METHODS 

In this section, we illustrate the simplified minimal 
model that we have selected to analyze the discrepan- 
cies between the LDA-f DMFT calculations on different 
(d-only vs. dp) basis-sets for the prediction of the orbital 
polarization and the Fermi surface properties of corre- 
lated multi-orbital systems. As discussed in the Introduc- 
tion, we know from preceding d-only and dp studies on 
bulk NickelatesfHI and Ni-based heterostructures [TtI [18| . 
that important discrepancies in the theoretical predic- 
tions emerged already when considering the simplest case 
of two correlated (namely, eg) orbitals 34] 

Hence, here, we set up two models with cu- 
bic/tetragonal symmetry to address precisely this issue. 
The models are obviously inspired by the realistic cal- 
culations of nickelate heterostructures, i.e. with a focus 
on the orbital polarization in a quasi two-dimensional 
geometry. However, as we are presenting DMFT calcula- 
tions, the "dimensionality" comes into play only via the 
fc-summations defining the local quantities, and, there- 



fore, the numerical results at the DMFT level will depend 
mainly on the kinetic energy of the two Cg orbitals, and 
on its relation w.r.t. the hybridization terms. Hence, the 
knowledge emerging from our analysis will also be used 
to clarify the general trends of the multiorbital physics 
at the DMFT level. 

Specifically, we will consider two models: 

1. The first one is the d-only two-orbital model (FiglT] 
upper panels) with the two orbitals representing 
two eg orbitals. 

2. The second model is a dp four-band model (FigH] 
lower panels) which, in addition to the Nickel 
3d — eg states incorporates two oxygen 2p states 
explicitly. 

Aiming at understanding the difference in orbital po- 
larization between the two models for a given electronic 
configuration, we have simplified each model to the min- 
imal set of parameters necessary to describe a physically 
realistic situation, without too many material specific de- 
tails. The two models are described below: 



d-only two-hand model 

For the eg only model in its simplest version we 
have only two parameters: Hopping amplitude tdd and 
initial crystal field potential /S.cf- The hopping was 
parametrized by means of cubic harmonics with the well 
known table of Slater and Koster [36] while the crystal 
field potential Acf accounts for the on-site energy dif- 
ference between the two eg states (i.e. it corresponds to a 
tetragonal distortion, a compression or elongation along 
one of the cubic axes, of the ligand octahedron). 

In addition, we pose a restriction to the filling, chosen 
to be one electron in the two effective eg states, corre- 
sponding to the downfolded states which are by no means 
of pure d-character. In fact, for the nickelates the two 
bands are the anti-bonding hybrid between the on-site-d 
and the ligand-p. The effective Hamiltonian in momen- 
tum space reads: 



■ cos ky) —^tddicos kx — cos ky) 
^tddicosk^ - cos ky) -Uddicosk^ + cos ky) + Acf 



|trfd(cos kx 

2 



(1) 



where the diagonal entries represent the hopping from a 
x^—y-^ (3z^— r^) to ax^— 2/^ (3z^—r^) in the next unit cell 
and the off diagonal elements represent the non-local Cg- 
Cg hybridization (cf. upper inset of FiglT|). Note in case 
of the heterostructure there is a quasi two-dimensional 
hopping geometry without hopping along the c-axis. This 



is why the 3z^ — has smaller hopping amplitudes and, 
consequently, a smaller bandwidth. A related two-band 
model Hamiltonian, but without hybridization between 
the orbitals, has been studied in DMFT in Ref. ^33], with 
reference to the cases of BaVSa and Na2,Co02. 

On the upper panels of Fig [T] we show the correspond- 
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FIG. 1. Band-structure and DOS of the two-band d-model (upper panel) and the four-band dp-model (lower panel) for two 
different values of the crystal field splitting Acf between the two d-orbitals, i.e. Acf = (left panels) and 0.5 (right panels). 
The orbital character is denoted by the following color-coding: black for the first d-orbital {x^ ~y^), red for the second (3z^ — r^), 
and yellow for the p-orbitals. In the inset, the hopping processes for the two models are visualized. 



ing band structure and single particle density of states 
(DOS). The color code denotes the orbital character: 
black and red represent the x'^ — y^- and — r^-orbital, 
respectively. While locally the two orbitals are eigen- 
states of the tetragonal point group and do not hy- 
bridize (pure black/red color of the DOS) the Cg orbitals 
obviously hybridize non-locally along certain directions 
(dark-red color, e.g. at the X point of the Brillouin 
zone). For the model we fixed the hopping amplitude 
to a value tdd — 0.6eV, which leads to a bandwidth com- 
parable to that of NMTO downfolding for the nickelate 
compounds [TtI. [T8| . The crystal field is the parameter 
for which we study several cases ranging from — 0.5eV to 
0.5eV. In the upper panels of Fig. [T] we show two cases 
^CF — {O.OeV, 0.5eV} as examples. The value of the 
local splitting can also be found in the difference of the 
center of mass of the DOS (drawn as dashed black/red 



lines). 

dp four-band model 

The second model we treat is in a way the "upfolded" 
version of the two-band model just described. We added 
two oxygen ligand orbitals in the quasi two dimensional 
geometry of the model which can be thought of as one 
and one Py orbital on each oxygen site. This is the 
minimal dp-model that can be constructed with a re- 
alistic dp-configuration in a quasi two dimensional cu- 
bic/tetragonal symmetry. For this upfolded model addi- 
tional parameters have to be assigned and chosen to be 
in correspondence with the two-band model: the param- 
eters (hopping amplitude and orbital splitting) of the d- 
only model should, as effective parameters, be derivable 
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from the parameters of the four-band model. Specifically, 
for the orbital splitting we remark that the formerly la- 
beled "crystal field" splitting is, in fact, a "ligand field" 
splitting which can be in principle decomposed into an 
electrostatic Madelung potential and a dp-hybridization 
splitting. Moreover, the dd-hopping processes are, in fact 
mediated by oxygen p-orbitals and we do not have a di- 
rect c?d-hopping in the four-band model. Keeping this in 



mind, we choose the parameters (dp-hopping amplitude 
and on-site d/p energies) in a way that the two bands at 
the Fermi energy (which are the anti-bonding dp-hybrid) 
reproduce the two bands of the previous model. We 
stress already at this point, however, that our subsequent 
results are robust with respect to the details of these 
parameters. The dp-Hamiltonian in momentum space 
reads: 



iff 



-i\/3tp 



;sm 



2tdd + A 



CF 



-it 



pd ' 



iVStpd sin(^) 



(Sin(^) 

En 



itr 







pd 



sm 



itpd sin 



(^) 







(2) 



where tdd and p are the parameters of the correspond- 
ing two-band model and tpd = 3 • tdd was chosen to yield 
the same bandwidth for the antibonding bands as in the 
previous model. The energy of the p-states was fixed to 
Sp = — 2.5eV similar to the position of the p-bands in the 
nickelate systems. 



DMFT for dp-models with projections 

Extended models including ligand p-states explicitly 
have a structure like our 4-band Hamiltonian Lo- 
cally, i.e. integrated over the Brillouin zone, such Hamil- 
tonians have the form: 



H\'Si{R = Q) = 



rrloc, 



H 



hyb. 

dp 



H 



hyb. 

dp 



rrloc 



(3) 



In the lower panels of Fig. [T] we show the band- 
structure of the four-band dp-model for the same crystal 
field splitting as for the two-band model. In the band- 
structure and DOS plots the additional yellow color en- 
codes the p-character. 



The local basis is typically chosen in a way that the 
Hdd and Hpp blocks are diagonal after k-integration so 
that the states can be labeled by a good local quantum 
number in the respective subspaces, such as the crystal 
field labels (see DOS plots in Fig. [1]). In such a basis, 
the local Coulomb {U) matrix of the interacting part of 
the Hamiltonian is then defined (in its SU(2)-invariant 
"Kanamori") form as 



h.C.) 



for the d-orbital sector. Here, U denotes the interac- 
tion parameter between two electrons in the same d- 
orbital, U' the interaction between electrons on different 
d-orbitals and J is the Hund's coupling; a, b index the two 
orbitals and a the spin. Note that the specific interaction 
values for the DMFT calculations (see next section for 
details) have been chosen in order to reproduce a typical 
correlated metallic physics situation. Obviously, due to 
the stronger localization of the d-orbitals in the dp mod- 
els, the corresponding values of the local interaction on 



the d orbitals have been correspondingly enhanced (we 
assumed here a factor two for the parameter U'). We re- 
call that -for the dp case- the multi-orbital Hubbard in- 
teraction could include, besides on-site d- and on-site p- 
interactions, also possible dp-interactions. This is, how- 
ever, not the topic of the present study. Further, for the 
dp model in LDA-I-DMFT we have to face the so-called 
problem of double counting correction (DC)0, IstI . [s^: 
This does no longer correspond to a simple total energy 
shift and, hence, cannot be "absorbed" in the chemical 
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potential. The DC for the dp-models corresponds to a 
renormalization of the energy difference between d and p 
states. For our models we have used the DC suggested 
by Anisimov [37| . given by: 

Udd ^ [U + U'{Nd - 1) + (C/' - J){Nd - l)]/{2Nd - 1) 

(4) 

where Nd is the number of c?— orbitals, and Udd is the 
average local interaction between these d-orbitals; n^^^ 
denotes the occupation of the d— orbitals as calculated 
from the k-resolved non-interacting (LDA) Hamiltonian 
discussed above. 

The self-consistent DMFT loop, which includes the so- 
lution of an Anderson impurity problem for the corre- 
lated subspace at each step is done as follows: 

1 . The first step is the calculation of the k-integrated 
Green function on the full dp basis set: 

Gtni^) f d?k [(c. + - iJ(k) - I]fuii(w)]-^ 

^BZ JbZ 

(5) 

where G, S, and H are matrices on the dp-basis. 

2. Next, we extract the dd-block of the local Green 
function: 

G',°;-(L.) = {G^°-H}U,_biock, (6) 

i.e. we project it onto the d-subspace. We stress 
here, that due to the dp hybridization encoded in 
the Hamiltonian -ff (k) and the inversion of Eq. [5] 
the information about the p-ligands is not lost but 
captured by G^^'^-{uj). 

3. Now, in complete analogy to DMFT for d-states 
only, we calculate the Weiss field for the impurity 
model (only on the d-subspace): 

where denotes the DMFT dd self-energy. 

4. With G°{ijJ) we solve the auxiliary impurity prob- 
lem (see below), obtain a new G|£f' and with the 
inverse of Eq. © a new T.^f^^{uj). 

Finally the self-consistent loop is closed by compar- 
ing both the new and the old 'S^^^'^{uj) and the 
new and the old d/dp-density and iterating until 
convergence. 

As impurity solver for our DMFT calculations, we 
have used the continuous time quantum Monte-Carlo 



(QMC) algorithm in the hybridization expansion (CT- 
HYB)[39i, i40i], which allows also for the treatment of the 
spin-flip and pair-hopping terms of our Kanamori Hamil- 
tonian (Eq. S]). In the CT-HYB the time evolution is 
calculated using the local interaction, which makes the 
local Hilbert space growing exponentially with the num- 
ber of orbitals. While for density-density interactions 
this problem can be mitigated along the line of ^3S. 4l|, 
for the more realistic Kanamori interaction a set of quan- 
tum numbers, called PS, leads to a more efficient algo- 
rithm [i^l . We have used this to perform all calculations 
presented in this paper with the full SU(2)-symmetric 
Hamiltonian. 



III. RESULTS: d VS dp CALCULATIONS AT 
"QUARTER FILLING" 

In this section, we compare our DMFT results for the 
d-only model at quarter-filling (ud = 1) with that of 
the corresponding dp model including also two 2p oxy- 
gen orbitals below the Fermi level (with a total density 
TT-tot — nd + Up — 5) . Specifically, in the next two sub- 
sections, we discuss how the final orbital polarization of 
the d-only (Sec. IIII A|) and dp calculations (Sec. IIII Bp 
evolves when varying the values of the initial crystal field 
splitting AcF between the d-orbitals. While the model 
cases studied here do not aim at a realistic description 
of a specific system, the results of this section can be 
qualitatively related, depending on the value of Acf = 
£^^2_j.2 — e.j.2_y2, to the physics of bulk Nickelates[T9| 
and of Ni-bases heterostructuresflTl fisl. Is^. In fact, (i) 
the nominal charge of these systems corresponds also to 
one electron in the outer two Ni-bands and (ii) the bulk 
Nickelates are typically characterized by negative values 
of Ac F , because of the tetragonal distortion along the z 
axis. In the Ni-based heterostructures instead, the local- 
ization effects in the z direction, as well as the epitaxial 
strain due to the substrate, induce positive values for 
Acf- 



A. DMFT results for the d-only model 

In this subsection, we analyze the results for our d- 
only two-orbital model at quarter filling (rid =^ 1), as a 
function of the initial crystal field splitting Aq'f- the 
corresponding results for the orbital occupation without 
interaction (which carry the label LDA, as this would 
correspond to the LDA in a realistic calculation) and 
with the interaction (computed with DMFT) are shown 
in Fig. [21 As mentioned in Sec. II, the interaction values 
have been chosen in consideration of typical values for 
transition metal oxides systems: for the d only model, we 
adopted a value oiU' = U-2J = 4 eV, with J = 0.5eV. 

We start by briefly commenting the set of non- 
interacting ("LDA") data shown in Fig. [2l They display 
a monotonous dependence on the initial crystal field: the 
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FIG. 2. Orbital occupation of the d-only model with U' — 
4eV, J = 0.5eV {U = 5eV), /3 = lOOeV"^ at quarter filling 
(rid = I) as a function of the initial crystal field splitting 
AcF- The DMFT data (solid symbols) are compared with the 
corresponding non-interacting results (empty symbols; coined 
LDA). The arrows indicate the effect of the interaction which 
is opposite for Acf > 0(< 0). 



orbital occupation of the — monotonously increases 
upon increasing values of Acf, whereas the occupation 
of the 3z^ — orbital decreases. We note, however, that 
as the hopping terms (and, hence, also the bandwidth) 
for the two orbital are not equal - a quite general situ- 
ation for realistic material calculations ~ the orbital oc- 
cupation curves are not symmetric w.r.t. to Acf- This 
also implies that the situation where the two orbitals are 
equally occupied (i.e. no orbital polarization) does not 
occur at Acf = but for the non-interacting case only 
for Acf < -0.5eV. 

We discuss now the effects of the Hubbard interac- 
tion on the orbital occupations, as described by our 
DMFT(CT-QMC) calculations. While the occupation 
curves remain obviously asymmetric also in presence of 
the interactions, from the data of Fig. [2]we note that, for 
each orbital, the deviations w.r.t. to the non-interacting 
values are strongly dependent on Acf- In fact, for all 
data considered, we observe that the sign of the change 
in the orbital occupation w.r.t. its LDA value appears 
closely connected with the sign of Acf: For Acf > 
we observe generally an enhancement (reduction) of the 



occupation of the first (second). 



-y^ (3z^-r^) orbital. 



while for Acf < the trend is the opposite. 

This represents, hence, a very robust behavior for the 
changes induced by the interaction to the orbital polar- 
ization P, which can be formally defined analogously to 
Refs. ^,111 



P = 



(8) 



FIG. 3. Difference of the real part of the DMFT self-energies 
for the two orbitals of the d-only model with U' = 4eV, 
J = 0.5eV {U = 5eV), 13 = lOOeV-^ at quarter filling {m = I) 
extrapolated to i:i;„ (solid symbols) and aj„ — s- +oo 
(empty symbol, Elartree contribution to the self-energy) as 
a function of the initial crystal field splitting Acf. The 
huge enhancement of the difference between the self-energy 
at ljJ„ — > 0, and the consequent huge energy shift of the first 
orbital, marks the onset of the Mott insulating phase (shaded 
region on the left). 



ferent LDA-I-DMFT results previously obtained for the 
shape of the Fermi surface via specific d-only calculations 
for bulk nickelatesfigj and Ni-based heterostructuresflj 
ITsj . respectively: The difference in the results simply re- 
flects the different sign of the initial crystal field splitting 
(Acf), as estimated by the ab-initio calculations, whose 
final size {Agff) is always significantly magnified by the 
electronic interaction. 

The last statement can be formalized more quanti- 
tatively through the analysis of the corresponding self- 
energies presented in Fig. [3j here we show the differ- 
ences between the real parts of the DMFT self-energy of 
the two orbitals, i.e. Rel]3^2„,.2 (iwn) — Rel]2.2_j,2 (ia;„), 
evaluated in the limit frequency ujn ^ and w„ — > +oo, 
respectively. We recall that in the latter limit only the 
Hartree contributions to the electronic self-energy re- 
mains. Hence, the difference between the self-energies 
can be also explicitly written in terms of the electronic 
density as: 



In fact, such a trend explains, on more general 
grounds, the physics underlying the qualitatively dif- 



Rel]3z2_r2 (oo) — ReS^2_j,2 (oo) — {U — 5 J) (9) 
{n^2^y2 - n3^2_,,2) = (C/ - 5 J) P. 

Here, the last equality only holds at quarter filling where 
n^2_y2 + 113^2 _^2 — 1. Such a dependence on the (final) 
electronic density appears evidently in the corresponding 
data of Fig. [31 where the sign change 4J] contribution to 
the DMFT self-energy is found precisely when the orbital 
occupations become equal, i.e. for the negative values of 
Acf ^ — 0.2eV for which P = in Fig. O By compar- 
ing this to the previous results and discussions, it should 
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be clear that the high-frequency values of the self-energy 
do not represent the "crucial" parameter. Instead, the 
trends in the orbital polarization and, therefore, the pre- 
dicted physics is controlled by the low-frequency behavior 
of the self-energy: Up to the point to which the systems 
still behaves as a (correlated) metal, i.e. when a Fermi 
surface still exists, the "effective" crystal field determines 
the splitting between the two orbitals in the presence of 
electronic correlations. It is given by 

Ae// = AcF + ReI]3,2_,2(0)-ReS],2_j,2(0), (10) 

i.e. the original crystal field is corrected by the difference 
of the real self-energies for w„ — >■ 0. This makes evident 
the interpretation of the second set of data shown in Fig. 
[3J the sign of Rel]322„^2 (0) — Rel]^2_y2(0) exactly fol- 
lows that of the original crystal field Aci?, and confirms, 
from the microscopic point of view, the quite robust pic- 
ture of interaction effects always magnifying the size of 
the original crystal field. Of course the magnification ef- 
fects will quantitatively depend on many factors. Quite 
naturally, they will be stronger when the system is more 
correlated. A dramatic enhancement is found when the 
Mott metal insulator transition is approached, i.e. when 
AcF approaches the shaded area in Fig. [3l Here the low 
energy physics correspond to an empty (broader) x'^ — 
orbital and a half-filled (narrower) — orbital. 

The coherent picture emerging from these 
LDA-I-DMFT, which would also provide a general and 
solid framework for understanding the previous d-only 
calculations for bulk nickelates and heterostructures, 
appears however to be contradicted when extending the 
basis-set of the LDA-I-DMFT algorithm to include also 
the most relevant p degrees of freedom, as it will be 
shown explicitly in the next subsection. 

B. DMFT results for the (corresponding) dp model 

We repeat in this subsection the same analysis done 
above for the d-only model, but now for the calcula- 
tions performed in the enlarged dp basis set, to study 
the differences between the two cases. Specifically, we 
consider here the second model in Sec. 2, including ad- 
ditionally, to the two d{eg) orbitals, also two p orbitals 
lying well below the Fermi level. The total occupations is 
ntot = rip+Ud = 5. Because of the stronger localization of 
the d-orbitals in the dp case, larger values for U, U' and J 
have been considered: For the sake of generality [l^ , we 
have simply doubled the value of U' and J w.r.t. the cal- 
culation discussed in the previous section, i.e. U' = 8eV, 
J — leV {U ^ U' — 2 J) which correspond, anyhow, to 
plausible values for correlated orbitals in transition metal 
oxides. 

As in the previous subsection, we consider first the 
plot (Fig. |4]) of the occupation of the two d-orbitals as a 
function of the initial crystal field splitting A.cf, start- 
ing from the non-interacting case. While the qualitative 
behavior as a function of Acf appears similar as before. 




-0.4 -0.2 0.2 0.4 

FIG. 4. Orbital occupation of the d orbitals as a function 
of the initial crystal field splitting Acf for the four-band dp 
model with U' = 8eV, J = l.OeV {U = lOeV), /3 = lOOeV-^ 
and the filling is ritot = 5. The DMFT data (solid symbols) 
are compared with the corresponding non-interacting results 
(empty symbols). 



an important difference should be noted: In spite of the 
relatively large separation (~ 2eV) among d and p band, 
due to the hybridization between the d and the p mani- 
folds, the total occupation of the d orbitals is now much 
larger than before {rid ~ 1-7 1.8). Quite remarkably, 
according to our DMFT results of Fig. 21 such an en- 
hanced occupation of the "correlated" d orbitals survives 
basically unchanged also when switching on the local in- 
teraction. 

This fact has an obvious impact on the final results for 
the orbital polarization P, as now the sign of its change 
w.r.t. the non-interacting case are no longer related to 
Acf- Fig- IH shows that independently on the sign Acf, 
one always observes a reduction of the value of P, i.e. 
a net enhancement (reduction) of the occupation of the 
3^2 _ j,2 ^^2 _ jy2-j Qi-]3j^ai driven by the electronic inter- 
action. 

This is evidently confirmed by the analysis of the corre- 
sponding self-energies: As discussed before, the most im- 
portant piece of information here is enclosed in the zero- 
frequency extrapolation of the real part of the DMFT 
self-energies. Fig. [5] shows the plot corresponding to Fig. 
[3] but now for the four-band dp-model: here the interac- 
tion correction always reduces the initial crystal field, in 
agreement with the systematic depletion of the x'^ — 
orbital (reduction of P) observed in the whole parameter 
range considered. 

The results are therefore in qualitative disagreement 
with those obtained within the d-only model, at least 
in the region of positive values of the initial crystal field 
Acf- It is worth noticing that for positive values of Acf, 
this trend is consistent with the results of the dp calcu- 
lations for the Ni-based heterostructures of Ref. [s^: 
There, it was shown, that even when starting from a 
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FIG. 5. Difference of the real part of the DMFT self-energies 
for the two d-orbitals of the four-band dp model as a func- 
tion of the initial crystal field splitting Acf with U' = 8eV, 
J = leV {U = lOeV), P = lOOeV"^ and the filling is ntot = 5. 
Shown is the low and high frequency asymptotes, i.e. the ex- 
trapolation to a;„ ^ (solid symbols) and ij„ +oo (empty 
symbol, Hartree contribution to the self-energy). 



reducing the value of J from l.OeV to J = 0.5eV, and 
J = O.OeV {U' = U -2J = 8eV is kept fixed, instead). 

The combined analysis of the orbital occupation and 
of the self-energy (at w„ — )■ 0) results shows that the dp 
results of Fig |4] and [5] change qualitatively already for 
the case J = 0.5eV: the overall trend of strong reduction 
of occupation of the first x"^ — orbital disappears for 
a large region of values of the initial crystal field Acf- 
In fact, the results already at J=0.5eV would lead to a 
physical situation which is qualitatively similar to that 
predicted by the d-only model. As one can expect, the 
change with respect to the previous dp results becomes 
even larger when setting J = O.OeV, as both trends of 
orbital occupation and effective crystal field becomes ex- 
actly opposite, with the orbital — y^ occupation always 
increased by the interaction irrespectively of the value of 
the initial Acf- The sensitivity of the final dp results 
on the Hund's coupling is very illustrative and shows the 
crucial role of the Hund's exchange in this situations pGj . 



IV. RESULTS: THE ROLE OF THE d-ORBITAL 
OCCUPATION 



relatively significant orbital polarization for the x'^ — 
orbital at the LDA level, the polarization was always 
strongly reduced by the interaction, which obviously has 
bad implications on the possibility of actually manipulat- 
ing the Fermi surfaces of the Ni-based heterostructures 
in the "desired" cuprate-like way [13, • 

As a matter of fact, this discrepancy between d-only 
and dp results has now been found here by including the 
SU(2) symmetry of the local interaction on the Cg bands 
in both DMFT calculations. This confirms the hypothe- 
sis of Ref. j3^] that such a discrepancy in the theoretical 
predictions does not originate from some different treat- 
ment of the interactions between the calculations of Refs. 
[13, [iSj and [32] , but rather from some intrinsic difference 
in the calculations. The most evident systematic differ- 
ence is the filling Ud of the two eg-orbitals, which, due 
to the dp hybridization is strongly increased w.r.t. the 
quarter filling level of the d-only model. While this is 
rather obvious at the LDA level, we note that the oc- 
cupation of the d-manifold does not change much even 
in presence of the interaction, for typical choices of the 
double-counting term for the DMFT (see Sec. 2, for de- 
tails, and also Ref. [111). Quite interestingly, the possible 
role of an enhanced d— orbitals occupations in dp calcu- 
lations has been also recently addressed, for the different 
problem of the occurrence of the MIT[31[. Here, besides 
the role of Ud, or more precisely, connected with it, we 
note that the final outcome of the results, and, in partic- 
ular, the occurrence of qualitative discrepancy between 
d-only and dp-results, is strongly infiuenced by the value 
of the Hund's exchange J constant. From the numerical 
point of view, this is exemplified by the results presented 
in Fig. [21 where the dp calculations have been performed 



From the marked discrepancy between dp and d cal- 
culations, discussed in the previous sections, a general 
question naturally arises: Under which conditions can 
one expect to obtain qualitatively similar DMFT results 
for d-only and dp-calculations? Let us assume here that 
we have a most accurate estimate [ESI of the local inter- 
actions and especially the Hund's exchange J. One of 
the most basic differences between the d-only and the dp 
calculations examined is the different value of the occupa- 
tion of the d orbitals. This difference is observed already 
at the level of the non-interacting "LDA" model and is, 
remarkably, not strongly affected by the inclusion of the 
interaction: For the d-only model, the LDA+DMFT cal- 
culations have been performed, as usual, at fixed filling, 
while the value of d-orbital occupation for the dp cal- 
culations presents rather moderate oscillations around a 
much larger value of Ud ~ 1.75, quite independently of 
the parameter set {U,J,Acf) considered. Since the fill- 
ing of the correlated orbitals is a crucial factor to drive 
the systems towards a Mott-Hubbard MIT, it is logical 
to suppose that it is also one of the pivotal parameters 
to be considered when comparing LDA-f-DMFT calcula- 
tions on different basis-sets. 

To make our statement more quantitative, we have 
performed several additional DMFT calculations for the 
same d-only model as in Sec. 2, but now varying sys- 
tematically Ud from the original quarter-filling level up 
to the much higher Ud ~ 1-75 found in the dp calcula- 
tions. The results for the orbital occupations and the 
self-energy at the Fermi level are reported in Fig. [7] in 
the upper and lower-row panels, respectively. By grad- 
ually increasing the total occupation of the d-orbitals, 
and focusing on the most interesting regime of Acf > 0, 
one observes that the parameter region where the crystal 
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FIG. 6. Left panel: Orbital occupation of the d orbitals of the four-band dp model as a function of the initial crystal field 
splitting AcF- We used U' = 8eV, j3 = lOOeV"^, ritot = 5 and two different values of the Hund's exchange J — 0.0, 0.5 eV. 
The DMFT data (solid symbols) are compared with the corresponding non-interacting results (empty symbols). Right panel: 
Corresponding data for the difference of the real part of the DMFT self-energies for the two d-orbitals of the four-band dp 
model extrapolated to tj„ — >■ as a function of the initial crystal field splitting Acf- The self-energy data are also compared 
with the corresponding DMFT data at J = 1.0, previously shown in Fig. 5. 



field "drives" the final result (i.e. where its magnitude 
is enhanced by the interaction), shrinks, being confined 
to higher and higher values of the initial A.cf- Looking 
at the self-energy plot for rid = 1-1, a qualitative change 
in the trend of the effective crystal field w.r.t. the non- 
interacting one is found for A^^^ 0.1. For Acf < ^cf 
the trend for the effective crystal field is opposite to that 
observed for quarter filling, i.e. the original crystal field 
Aqp is reduced by the interaction. A similar situation, 
though for a shghtly larger threshold (A^^;, ^ 0.15), is 
also observed for the changes in the orbital occupations, 
w.r.t. the non-interacting ones. By further increasing 
to 1.25 (second panels of Fig. [71)), the threshold A^^ 
is already shifted at the border of the parameter region 
considered (AJ.^ ^ 0.5). Finally, if one performs the 
DMFT c?-only calculations for a similar = 1.75 as in 
dp-model, one indeed finds similar results as for the dp- 
model: Both, the effective crystal field and the changes 
in the relative orbital occupations in Fig. [71 = 1.75 
resembles those of the dp-model in Figs. [31and[5l 

Let us stress that due to the different bandwidths of 
the two eg orbitals, no symmetric behavior between the 
regions of positive and negative Aqf can be expected. 
In fact, in the region Aqf < (i.e. the less relevant one 
in the perspective of the Ni-based heterostructures), the 
crystal field "enhancement" at quarter filling was much 
larger than the corresponding one for Acf > (see Fig. 
121), also as an effect of the closer proximity of the MIT 
in this parameter region. Hence, while the effects of an 
increased Ud are always at work, here they are just visible 
as a mere weakening of the quarter fillings trends rather 
than as of an inversion of them: The gray shadow area, 
marking the onset of the MIT for one of the two orbitals 
i.e., the (almost) half- filled one, is shrinking in the Ud = 



1.1 plot, and essentially disappears at Ud = 1.25. At 
higher densities, one also observes that the large size of 
the (negative) effective crystal field for Acf = —0.5 is 
visibly reduced w.r.t. the case of quarter filling. 

These results demonstrate the important role played 
by the density of the d-orbital manifold for determining 
the final DMFT results on different basis sets, at least 
for the important aspect of the Fermi surface properties. 
The starting value of nd will decide, whether, given a 
correct ab-initio estimate of the interaction parameters 
of the multi-orbital Hubbard Hamiltonian, the physics 
of the interacting system will be driven by the (original) 
crystal field, or rather by the Hund's rule tendency for 
equally occupied orbitals. Therefore, it will be a priori 
quite hard to reconcile DMFT calculations performed on 
different basis-sets, without considering the correspond- 
ing occupation of the correlated orbitals manifolds. 

From the perspective of actual material calculations, 
the strong dependence of the final LDA+DMFT results 
for the Ni-based heterostructures is now much easier to 
understand. In fact, the quarter-filling physics particu- 
larly favors the crystal field-dominated physics, since the 
Hund's exchange is weaker in a system with one electron 
(on average). Hence, if the dp hybridization can drive 
the system away from this regime (as is the dp model 
considered here with an average filling of Ud ^ 1-75), the 
Hund's exchange will eventually prevail over the crys- 
tal field effects. The same interpretation can likely ex- 
plain, why the results of LDA+DMFT calculations for 
the Fe-based superconductors do not display such cru- 
cial dependence on the basis-set considered. There, the 
five orbitals of the Fe-Sd manifold are characterized by 
a very small value of Acf ~ 0.2eV (in comparison to 
the typical values of leV observed in many transition 
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FIG. 7. Upper row: Orbital occupations of the two orbitals of the d-only model with U' = 4eV, J — 0.5eV {U = 5eV), 
/3 = lOOeV"^ for different fiUing na = 1.1, 1.25, 1.75 as a function of the initial crystal field splitting Acf- The DMFT data 
(solid symbols) are compared with the corresponding non-interacting results (empty symbols). Note that a filling of = 1.75 
roughly corresponds to the filling of the two d-orbitals in the dp models considered in the previous sections. Lower row: 
Difference of the real part of the DMFT self-energies for the two orbitals extrapolated at aj„ — (solid symbols) as a function 
of the initial crystal field splitting Acf for the corresponding set of data. 



metal oxides). This corresponds to a situation of five 
partially filled (n^ ~ 6) correlated orbitals all very close 
energetically at the Fermi level, i.e. the most favorable 
playground for strong Hund's exchange processes. Hence, 
this case is well inside in one of the two regimes, and pos- 
sible moderate differences in the initial orbital occupation 
in the c?-manifold in different basis-set will not result in 
a qualitative discrepancy between different LDA-I-DMFT 
calculations. 



V. CONCLUSIONS 

We have thoroughly compared two models for transi- 
tion metal oxides: a c?-only model with two effective d- 
orbitals and a dp-model with two d- and two p-orbitals. 
On the single-particle level without Coulomb interaction, 
the two models show the same low-energy physics and 
band-structure. However if the Coulomb interaction is 
taken into account by means of LDA-I-DMFT, this is not 
the case any more. 

The main reason for this discrepancy is the number of 



d-electrons. In the d-only model there is on average one 
d-electron per site. For the dp-model on the other hand, 
the bands are filled with altogether 5 electrons per site: 
Without hybridization the two p-orbitals at lower energy 
would be completely filled, and the d-orbitals in the vicin- 
ity of the Fermi level would have the same filling (one 
electron per site) as for the d-only model. The d-orbitals 
hybridize however with the p-orbitals, and, hence, there is 
some admixture between the orbitals leading altogether 
to a larger filling (^ 1.7) of the d-orbitals for the dp- 
model. 

This different filling has dramatic consequence for the 
correlated solution of the two models. For the d-only 
model the Coulomb interaction enhances the initial crys- 
tal field splitting, leading to a situation with one d- 
electron in the lowest-lying d-orbital. Depending on the 
original crystal field (e.g., if Acf > 0), this can be the 
— y^-orbital, which hence dominates the low energy 
physics and Fermi surface topology. This kind of physics 
can result in a very similar occupation of the x'^ — y^- 
orbital and Fermi surface in Ni-based heterostructures 
as in cuprates. 
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In contrast, for the rfp- model the Hund's exchange fa- 
vors a situation with a more even occupation of the two 
orbitals, an enhancement of the local magnetic moment. 
In this situation, the initial crystal field splitting is actu- 
ally reduced. Correspondingly, the initial disproportion- 
ation in the occupation of the two orbitals is shrinking. 
We have verified that the d-only model shows the same 
kind of Hund's physics if we artificially consider the same 
d-electron occupation as in the dp-model. 

Our findings have important consequences for 
LDA-f DMFT calculations of transition metal oxides and 
related correlated materials: The physical results for d- 
only LDA-I-DMFT calculations can be dramatically dif- 
ferent from calculations which also include the oxygen p- 
orbitals. Since the additional p-orbitals lead to a different 
d-fiUing the local interaction can result in the aforemen- 
tioned dramatically different local DMFT physics. Even- 
tually, experiments (such as, e.g., ARPES or X-ray ab- 
sorption spectroscopy or orbital reflectometry) will show 
which of the two physical scenarios is correct. 

Furthermore, as the dp-calculation is generally far 
away from an integer filling of the interacting d-orbitals, 
the observation of Mott- Hubbard metal-insulator transi- 
tions becomes more difficult [30,, i31|]. While one can ex- 



pect non-local correlations beyond DMFT to be of impor- 
tance for two-dimensional Mott insulators, not showing 
a Mott insulating phase in specific other cases might be 
a deficit of the dp-calculations - or, at least, of the way 
such dp-calculations are performed nowadays. 

Let us emphasize that in other situations the physics 
of the d- and dp-calculation can be much more similar. 
This was e.g. observed for iron pnictides. In these mate- 
rials, we have a large number of d electrons for both dp- 
and d-only calculation so that in both cases the physics 
is dominated by Hund's rule forming a local magnetic 
moment. 
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